Local topological and chiral properties of QCD 
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To elucidate the role played by instantons in chiral symmetry breaking, we explore their properties, in full QCD, 
around the critical temperature. We study in particular spatial correlations between low-lying Dirac eigenmodes 
and instantons. Our measurements are compared with the predictions of instanton-based models. 



We have examined the local topological struc- 
ture of full QCD and its impact on the physics 
of chiral symmetry breaking and restoration, by 
comparing the local topological structure ob- 
tained by improved cooling yj], with the lowest 
eigenmodes of the Dirac operator on the same 
(uncooled) lattices. We use a set of dynamical 
finite temperature N t = 12 (MILC@]) lattices 
spanning the chiral phase transition, and have ex- 
tracted the lowest 8 eigenmodes and eigenvalues 
of the staggered Dirac operator. 

To summarize our findings, we see: 

• agreement with the Banks-Casher relation- 
ship using the density of eigenvalues near zero 
(not discussed in this writeup). 

• good correlation (~ 70% level) between the 
spatial structure of zero modes and instantons. 

• some space-time asymmetry in the local topo- 
logical susceptibility above T c (and not below), 
the underlying mechanism of which is under fur- 
ther investigation. 

1. LOCAL STRUCTURE 

To illustrate the relationship between instan- 
tons and the zero modes of the Dirac operator we 
show in Figures 2 and 3, the topological charge 
density FF(x), on a configuration obtained after 



150 sweeps of improved cooling. The timeslice 
shown happens to contain part of an instanton 
and anti-instanton, and we plot isosurfaces of pos- 
itive and negative values of FF(x). 

On this cooled configuration, we identify the 
lowest eigenmode of the Dirac operator; an isosur- 
face of the magnitude of the eigenvector is plotted 
along with FF(x), and is shown in Figure 2. We 
see that on this smooth configuration in which 
the UV fluctuations have been removed, \ipip(x)\ 
follows FF(x) exactly, showing that on contin- 
uum like configurations, the zero mode "tracks" 
the topology, as expected from continuum argu- 
ments. 

Next, Figure 3b compares the uncooled zero 
mode to the cooled topological charge; we see sur- 
prisingly good correlation, even after many (150) 
cooling sweeps. Cooling identifies the dominant 
instanton — anti-instanton (I- A) pairs and the en- 
semble correlation between the uncooled Dirac 
mode and cooled topological charge density is 
about 70%, after 150 sweeps of improved cool- 
ing, on configurations containing one or more TA 
pairs. This validates a posteriori the improved 
cooling process; the topology seen by the Dirac 
zero mode is largely the same as that seen by the 
cooled topological charge. 
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In Figure 3a, we show isosurfaces of ^^(x), 
which takes value ±1 for a right- or left-handed 
eigenmode respectively. Yellow and blue surfaces 
indicate right and left-handedness. We see what 
amounts to the chirality flip as quarks "tunnel" 
from an instanton to an anti-instanton. 

2. ABOVE T c 

At T ~ 1.25 T c , all configurations have Q=0, as 
X — > sharply. It is very difficult to say that there 
are instantons and anti-instantons. After a min- 
imal degree of smoothing (~20 cooling sweeps) a 
fit to non-interacting (anti)instantons or calorons 
does not converge. Before that, the configura- 
tions are too rough to attempt any topological 
identification, and we cannot really say for sure 
that there are instantons at all. 

We can however try to characterize the fluctu- 
ations of FF(x). In instanton liquid models j3), 
instanton — anti-instanton dipoles are thought to 
orient in the timelike direction. Inspired by this 
theoretical picture we examined our data in the 
following way||]. 

We slice the lattice up into smaller sublat- 
tices, taking slabs of a given thickness d. While 
the slabs have a given thickness in direction //, 
they span the lattice in the other three direc- 
tions, fi is chosen either along the time axis or 
one of the spatial axes, so that our slabs span ei- 
ther space or time. Within the slab we compute 
J[FF(x)] 2 d 4 x = Qg lab , then average over all pos- 
sible slabs of thickness d which are of the same 
orientation w.r.t. space or time. 

If instantons and anti-instantons are randomly 
placed in an isotropic lattice, we expect that, 
up to a volume normalization factor, spatial and 
temporally oriented slabs of the same thickness 
d should have the same value of (Q 2 )d- This is 
what we see on the lattices below T c . 

If dipole pairs have formed in the timelike direc- 
tion, we would expect (Q 2 )d to plateau for large 
d for slabs which span the time-like direction (fi 
in the spatial direction). This is because we are 
always adding a Q = pair to the slab as we in- 
crease its thickness, and this does not increase Q 2 
in the slab. For spatially spanning slabs on the 
other hand, the addition of charge is more or less 



random. 

What we see above T c , after a small amount of 
cooling (less than 20 sweeps), is the behavior just 
described: a plateau in {Q 2 )d for d=spacelike ori- 
ented slabs (time-spanning), and no such plateau 
in the corresponding G?=timelike (space-spanning) 
slabs. Although it is quickly washed out by cool- 
ing, the asymmetry is quite clear and unmistak- 
able. What is less clear is its interpretation. 

0.6 I 1 1 1 1 1 1 1 1 1 1 1 1 1 

£ = 5.85 r [ " 



0.4 — 

A 



O" 

V 



0.2 — 



0.0 




5 10 

Slice volume 

Figure 1. (Q 2 ) versus volume: The lattice is 
sliced into sections of given volume, and (Q 2 ) is 
computed on these. The solid curve has time- 
spanning slices (d in the spatial direction) and 
the dashed curve has space-spanning slices (d is 
in the timelike direction). If there are dipoles in 
a particular direction (timelike it would appear), 
then (Q 2 ) should plateau when slices are taken so 
that they span that direction. Notes: the spatial 
axes have been rescaled to match the temporal 
dimension (N t = 12), and Q 2 is measured after 5 
cooling sweeps. Below T c , both curves are identi- 
cal within statistics. 



It is important to note that we see this asym- 
metry for quenched SU(2) as well, where it is not 
predicted to be by instanton liquid models. 

Since obtaining this result we have generated 
synthetic configurations with a controlled topo- 
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logical charge density, and computed (Q 2 )d on 
these. On synthetic configurations we see: 

• For randomly placed point charges, q(x) = 
±<5(x — Xi), there is no asymmetry. This cor- 
responds to the data below T c . 

• With randomly oriented dipoles having point 
charges, we see a trivial asymmetry of geometric 
origin. There is a plateau in both space and time 
slices, and the ratio of the plateau values is equal 
to the ratio of the slab sections (i.e. 2). 

• For randomly placed calorons, we see a re- 
verse asymmetry which shows dependence on the 
radius of the caloron. By reverse asymmetry, we 
mean that it is the space-spanning (d timelike) 
slabs which plateau at a lower value, the reverse 
of what is seen in Figure 1. 




Figure 2. Isosurfaces of topological charge den- 
sity FF(x) (red[+] and blue[-]) and tpo(x)i(>o(x) 
(yellow) using the lowest eigenmode ipo- Here 
?po(x)?po(x) is computed on the cooled configu- 
ration. 
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Figure 3. Isosurfaces of topological charge den- 
sity FF(x) (red[+] and blue[-]) and tpo(x)'ilJo(x) 
(yellow) using the lowest eigenmode Vo- Here 
ipo(x)ipo(x) is computed on the original uncooled 
configuration. Above (a), chiral components; Be- 
low (b) | tp ip (x)\ 



